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Recently a new approach to calculate the nuclear potential from lattice QCD has been proposed. 
In the approach the nuclear potential is constructed from Bethe-Salpeter (BS) wave functons 
through the Schroedinger equation. The procedure leads to non-local but energy independent 
potential, which can be expanded in terms of local functions. In several recent applications of 
this method, local potentials, which correspond to the leading order (LO) terms of the expansion, 
are calculated from the BS wave function at £ ~ MeV, where E is the center of mass energy. 
It is therefore important to check the validity of the LO approximation obtained at £ ~ 0. In 
this report, in order to check how well the LO approximation for the NN potentials works, we 
compare the LO potentials determined from the BS wave function at £ ~ 45 MeV with those 
at £ ~ MeV in quenched QCD. We find that the difference of the LO potentials between two 
energies are not found wihin the statistical errors. This shows that the LO approximation for the 
potential is valid at low energies to describe the NN interactions. 
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1. Introduction 

The nucleon-nucleon (NN) potential is widely used in nuclear physics. Once the NN potential 
is known, one can, in principle, determine various structure of nuclei, by simply solving the cor- 
responding Schrodinger equation. In the past few decades, several different NN potentials, such 
as phenomenological potentials determined from the fit of NN scattering data with ;^^/dof ~ 1 
at Tiab < 300 MeV [|l|, |2[ ^ or potentials based on effective field theory [Q], have been proposed. 
Although these potentials can reproduce the scattering data to a quite good accuracy, they have a 
drawback that they need a large number of parameters to describe the phase shift of the NN scat- 
tering (AVi8 needs 40 and ChPT in N^LO needs 24 parameters. ). In this respect, the NN potential 
from Lattice QCD proposed in Ref. [jsj, |6l, 0] has an advantage: It requires only a few fundamental 
parameters of QCD, the gauge coupling constant and quark masses mu,md,ms,- ■ ■, so that the 
method can be also applied to hyperon systems {NA,NL . . . ) [|8|], for which only a limited number 
of experimental information is obtained so far. 

Recently Hadrons to Atomic nuclei from Lattice (HAL) QCD collaboration is formed to study 
various aspect of baryon-baryon potentials based on the first principle of QCD. In the approach 
by the HAL QCD collaboration, the non-local but energy-independent potential U{r,r') is first 
constructed from the Bethe-Salpeter (BS) wave function via the Schrodinger equation [^, ^, ^. 
The non-local potential U{r,r') can be described in terms of local functions using the derivative 
expansion, whose leading order contribution gives the local potential. In several applications of 
this method, the leading order local potentials have been evaluated at zero energy in the center 
of mass frame. Therefore, in practice, it is important to know at which energy range the leading 
order local potential is accurate enough to approximate the non-local potential U{r,r'), which is 
faithful to the scattering data by construction. In this report, we extract the leading order local 
potentials at non-zero energy from the quenched lattice QCD simulation, and compare them with 
those previously obtained at zero energy. A difference between them gives an estimate of higher 
order correction in the derivative expansion. 

This report is organized as follows. In section |2|, we give a brief review of the method to 
extract the NN potential in Lattice QCD using the derivative expansion. In section ^ we compare 
the leading local potentials between zero and non-zero energies. We have found that a difference 
between them is small compared to statistical errors. In section ^ we consider contaminations to 
the potentials from excited states, which become manifest at large separation where the potential is 
expected to vanish. Section ^ is devoted to summary and conclusion. 

2. NN potential from Lattice QCD 

The non-local potential U {f,r') is constructed from the equal-time Bethe-Salpeter (BS) wave 
function (x, k) through the Shrodinger equation ^ as 

(A+A:^) ^{r;k) =mN I r' U{r,r) ^(f ,k), (2.1) 

where "k" denotes the "asymptotic momentum", which is related to the total relativistic energy 
W asW = iJ rn^ + k'^. The derivative expansion up to V, together with various constraints from 
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symmetries leads to the conventional form of the NN potential at low energies widely used in 
nuclear physics : 



U{r,r')= Vo'(r)+Vi(r) (ai • 02) + 5i2 + V/^W L ■ 5+ ^(V2)' 



(2.2) 



where r = |r|, Sn = 3(ai ■r){d2-r)/r^ — 01-02 is the tensor operator, S = {01 + 02) /2 is the total 
spin, L = r X p is the orbital angular momentum, and / = 0, 1 is the total isospin. Note that L ■ S, 
which is ff{V), is of next-to-leading order in the expansion. Eq. (2A) with (2^) successively 
determines local functions V((x) (A = 0, o, T,LS, ■ ■ •)■ 

The BS wave function on the lattice with the lattice size L is defined by 

<^'{r,k) = ^Y^P'„p{0\np{r + x)pa{x)\B = 2-,W), W = 2^ml + k\ (2.3) 

X 

where W is the total energy of the two nucleon system in the center of mass system, and P^^ denotes 
a projection operator to the spin singlet state {Pf^^ = (o'2)a./5) or triplet state {P^^ = (<7i)a,j8) for 
spinor indices a and j8 . Local composite operators for the proton and the neutron h and p are given 
by 

n^{y) = £abc {uaiy)CY5dh{y)) dcpiy), Paix) = Eabc {ua{x)Cysdb{x)) Uca{x), (2.4) 

where, a,b,c denote color indices, and C is the charge conjugation matrix. 

In this report, we consider potentials for the '^o state and the ^Si - state. In the case of 
state, the Schrodinger equation at leading order becomes 

(A + ;t2)0'^o(r;A:) = 2jUy'^"(r)0''^«(r;A:) (2.5) 

for the spin singlet channel {S = 0) with the reduced mass }X = 111^/2, where the wave function 
^ ^°{r,k) for the state is given by the projection P as 

0'*(r;^) = P<^\r,k) ^ £ <^\R[n\k). (2.6) 

■^^ ReO 

Here the summation over /? € O is taken for the cubic transformation group to project out the 
state. Then the central potential is easily obtained as 

where £(= is an effective kinetic energy in the center of mass system. 

For the spin triplet channel {S = 1), the Schrodinger equation at leading order becomes more 
complicated due to the mixing between and components by the tensor potential: 

{A + k^)^\r;k) = 2pL{V,{r) + VT{r)SnW{r-k) (2.8) 

where Vc(r) = Vq^C?") + V^(r). By two projections P and Q = \—P, the above Schrodinger equation 
is decomposed as 

P<^\r;k) PSi2(p\7;k) \ ( K(r) ^ _ A + ^ ^ ^^^^ 



\Q(t>Knk) QSn(l>Hnk) J yVrir) J 2^1 \Q<i>\r,k) 
from which Vdr) and VV(r) can be separately obtained[[T^, ^]. 
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3. Numerical Simulations and results 

3.1 Lattice QCD setup 

We employ the standard plaquette gauge action on a 32^ x 48 lattice with the j8 = 5.7 for 
quenched gauge configurations. Quark propagators are calculated by the Wilson quark action at 
K = 0.1665. This setup leads to the lattice spacing a^^ = 1.44(2) GeV (a ~ 0.137 fm) from nip, 



the spatial extension L = 32a ~ 4.4 fm, mj^ ~ 0.53 GeV and ~ 1.33 GeV [|ll[]. Quenched 
gauge configurations are generated by the heatbath algorithm with overrelaxation. Potentials are 
measured on configurations separated by 200 sweep. 4000 configurations are accumulated to obtain 
results in this report. These calculations are performed on Blue Gene/L at KEK. 

The BS wave function is obtained from the four-point correlator of nucleon operators in the 
large t region, 

G'^^\x,y,t,to) = (0|«/3(y,0^a(x,0/p«(?o)|0) =£A„(0|«/3(3;)^„(x)|B = 2;W„)e-^"('-''') 

n 

r^Ao{0\np{y)pa{x)\B = ^-,Wo)e-''''-'-"\ Wn=2^ml + kl (3.1) 
Here the source located att = tQ, ^pn{k)), is defined by 

/pn = P^'P'Pa' ito)Npi (to) , Pa = Ea.b^c (UaCYsDh) Oca , Np = £a,b,c {UaCysDh) D^p , 

where 

U{t) =Y^u{t,x)f{x), D{t)=Y^d{t,x)f{x), (3.2) 

X X 

where f{x) is source function, as will be seen later. By examining the t dependence of potentials, 
we see that ground state saturations for potentials are achieved att — tQ = 9. 

3.2 Periodic and anti periodic boundary conditions 

The periodic boundary condition (PBC) is imposed on the quark fields along the spatial direc- 
tions to obtain the NN potential at £ ~ MeV, while the anti-periodic boundary condition (APBC) 
is employed for the NN potential at £" 3 x (n/L)'^ /m^, which corresponds to 45 MeV in the cen- 
ter of mass system. For the PBC, we employ the wall source, i.e., f{x) = 1 , which enhance ground 
state of PBC, p = (0,0,0)7r/L. On the other hand, for the APBC, we employ four types of momen- 
tum wall sources, f{x) = cos{{+x+y + z)7t/L), cos{{—x+y + z)7t/L), cos{{+x — y + z)7r/L), and 
cos{{—x—y+z)7i/L), where these sources enhance the ground state of APBC, i.e., p = (1, 1, l)7r/L 
state. Here, we have imposed positive parity to the system by using the wall source with a cosine 
type instead of an exponential type. After the summation over the results of four sources, the 
representation is obtained. To improve the statistics for the PBC case, we locate four sources on 
different time slices on each configuration. 

To evaluate the central potential by eq. (2.7), we need to determine the first term, /IpL, 



which, in principle, is obtained through the relation Wo = 2y m'^ + k'^. Within statistical and sys- 
tematic errors, however, the values of turn out to be close to the free values. We therefore adopt 
the free values, E = k^/triN = MeV (PBC) and E = k^/rriN = 3{n/Lf/mN = 45 MeV (APBC), 
to calculate potentials in this report. 
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Figure 1: Comparison of local potentials with the PBC (blue open circle) and the APBC (red closed circle) 
at f — fo = 9 and r < 2.0 fm. The left figure shows the '5o central potential, while the center and the right 
represent the — ^ Di central and tensor potentials, respectively. 



3.3 Comparison of local potentials between two energies 

In Figure [l], we compare potentials at the leading order of the derivative expansion obtained at 
£■ ~ 45 MeV (red circles) with those at £" ~ MeV (blue circles), for the v'^°{r) (the left), K (r) 
(the center) and Vrir) (the right). All data are taken at f — ?o = 9, where grand state saturations 
for the potentials are achieved. From these figures we observe that the agreement of potentials 
between two energies is quite good for all cases within the statistical errors. We therefore conclude 
that the leading order contribution in the derivative expansion gives a very accurate approximation 
for the energy-independent non-local potential U{r,r') in the energy range from E = MeV to 
E = 45 MeV, in the case of the quenched approximation at nijc ^ 530 MeV. Note that the NLO 
contribution is absent for the spin-singlet channel while the NLO Vis potential exists for the spin- 
triplet. The agreement of Vc and Vj between two energies suggests that Visir) is sufficiently small 
below £" = 45 MeV, at least for nij^ ~ 530 MeV. 



4. Contamination from excited states at large distances 

Results in the previous section show that the local potentials obtained from the BS wave func- 
tion at £■ ~ 45 MeV agree well with those at £" ~ MeV. One may notice, however, that potentials 
obtained at £ ~ 45 MeV (APBC) deviate from zero at large r, where potentials are expected to 
vanish. These deviations are not statistical fluctuations , as seen in Fig.^, where the local potentials 
obtained with the APBC are plotted as a function of r at ? — fo = 3, 6, 9: Deviations of the potentials 
from zero at large distances have characteristic structures, which are most clearly seen at ? — ?o = 3, 
and the deviations decrease as f — to increases. These observations suggests that the deviations of 
the potential from zero are caused by contaminations of excited states. 

Here we discuss how contaminations of excited states to the BS wave function affects the form 
of the corresponding potential. Let us assume that the BS wave function extracted from the 4 point 
function is dominated by the grand state with a small contamination of one excited state as follows. 

^{r,t) = Y{7,ko)e-''" + Y{7M)e-'^'', {ko<h), (4.1) 

where W„ = lyjmj^ + kl, and i/A(r,/:o) is the BS wave function of the grand state with Eo = kQ/mj^/ 
while Y{r,ki) is the wave function of the excited state with E\ = k\/mj^. At sufficiently large r. 
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Figure 2: These figures are ^So central potential with the APBC at f — fo = 3 (left), 6 (center) and 9 (right). 
These potential deviate from zero at r > 2 fm. Deviations of potentials are decrease as f — fo increases. 



we assume that both wave function satisfy 



1 



(4.2) 
(4.3) 



By inserting eq. ( |4. l] ) into eq. ( |2.7[ ), we arrive at the following expression: 

1 A^{f,t) kl r^oo 1 ,2 ,2\ 



+ ■ 



where P{r) is given by 



P(f) 



IflN 



{k\-kl) P{r), 



(4.4) 



Y{r,ko) + \l/ir,h)e-"-^'-^°^ 



-t{W,-Wo) 



(4.5) 



We see that, even in the non-interacting region where the true potential V{f) vanishes, eq ( p^ 
leads to a non-vanishing "potential" as 



niN 



(4.6) 



The agreement of local potentials between two energies indicates that the effect of these con- 
taminations to potentials is smaller than statistical errors at short distance (r < 2.0 fm). Therefore 
the conclusion in the previous section remains true. Since the true potentials vanish at long dis- 
tance, however, the small effect due to the contaminations can become significant. Note that the 
APBC implies that not only the numerator but also the denominator of P(r) vanish at boundaries, 
so that P{f) could become large near boundaries. The consideration so far suggests that the devi- 
ations of the potentials from zero with the APBC at large r is caused by the contaminations from 
the excited states. Our choices of momentum wall sources creates not only the grand state with 
p = £(ibl,zbl,ibl) but also the excited state with p = J(ib3,±3,±3). Moreover the energy dif- 
ference in lattice unit is not so large: Wi — Wq ~ 0.23/a (~ 360 MeV) in this case. Therefore it 
is likely that the contamination comes from the (3,3,3) excited state. This observation, however, 
needs to be confirmed, which is now underway. 
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5. Summary and conclusion 

We have examined how well the leading order contributions in the derivative expansion of 
the non-local potential U (r,r') describe the NN interactions in the wide range of energy. We have 
compared the local NN potentials for the ^Sq state (the central potential) and for the ^S\-^D\ state 
(the central and the tensor potentials) obtained at ~ 45 MeV with those obtained at £ ~ MeV 
in quenched QCD at rrij^ ~ 530 MeV. We have found that differences of these potentials between 
two energies are very small. From this result we conclude that the leading order local potentials in 
the derivative expansion are good approximations for the NN potentials at energy up to 45 MeV in 
quenched QCD, and that the local potential constructed at £" ~ MeV can be used to investigate 
properties of the NN interaction at low energy. 

In the future it is important to apply the analysis in this report to the NN potentials in full 
QCD at lighter pion mass and to more general potentials including hyperons, in order to confirm 
the validity of the leading order local potential approximation for these cases. 
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